November 23, 2004 


Dr. Philippe Crane 
Origins Theme Scientist 
NASA Headquarters Code SZ 
300 E St S.W 
Washington, D.C. 20546 
Tel: (202) 358 - 0377 
Fax: (202) 358 - 3096 
Email: pcrane@hq.nasa.gov 


Final Progress Report for NAG5-9581 

Dear Dr. Crane, 

This letter is the final progress report for the SARA grant NAG5-9581 entitled “Improved 
Simulations of Astrophysical Plasmas: Computation of New Atomic Data” which 
began on 1 March 2001 and ended on 28 February 2004. 


1 Introduction 


Our research program is designed to carry out state-of-the-art atomic physics calculations 
crucial to advancing our understanding of fundamental astrophysical problems. We redress 
the present inadequacies in the atomic data base along two important areas: dielectronic 
recombination and inner-shell photoionization and multiple electron ejection/ Auger fluores- 
cence therefrom. All of these data are disseminated to the astrophysical community in the 
proper format for implementation in spectral simulation codes. 


2 Review of Progress 

2.1 Brief Summary 

The following is an initial, brief summary of the atomic physics research we have performed 
over the grant period, which was funded by NASA Space Astrophysics Research and Analysis 
Program grant NAG5-10448: 

• Benchmarking of our MCBP and R-matrix DR calculations with experimental mea- 
surements performed at the Test Storage Ring (5 iron ions, 4 papers); 
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• Systematic calculations and fitting formulas for DR of C-like, O-like, F-like, and Ne- 
like isoelectronic sequences (95 ions, 3 published papers and one more in preparation, 
tables of fitting coefficients on our web page); 

• Assesment of the exisiting fluorescence/ Auger database (1 paper); 

• Peripheral collisional (5 papers), photoionization (2 papers), photodetachment (3 pa- 
pers), and DR perturbative vs. R-matrix comparison (1 paper) studies. 

A more detailed description of each of these accomplishments is given below. 


2.2 General Considerations 

Our research program is designed to carry out state-of-the-art atomic physics calculations 
crucial to advancing our understanding of fundamental astrophysical problems. We have 
redressed the present inadequacies in the atomic data base along two important areas: di- 
electronic recombination and multiple fluorescence/ Auger decay of inner-shell vacancy states. 
All of these data have been disseminated to the astrophysical community in the proper format 
for implementation in spectral simulation codes. 

With the available funding, we were able to support a postdoctoral fellow for only two of 
the three years of the grant period. On November 1 st , 2001, we hired Dr. Oleg Zatsarinny, 
an expert in atomic collision calculations, and the former director of theoretical physics at 
Uzhgorod University, Ukraine. 

Dr. Zatsarinny thus worked at WMU for two years. During that time, he diligently performed 
DR calculations for 95 different ionic systems , as is detailed below. 

Besides this large comprehensive effort of computing along isoelectronic sequences, our other 
two main areas of research have been 1) the benchmarking of theoretical and experimental 
DR rate coefficients for several iron ions, and 2) the assessment of the fluorescence/ Auger 
data base in current use. Further work in collisional excitation, photoionization, and pho- 
todetachment made up smaller peripheral projects. 


2.3 DR of O-like Fe XIX: Laboratory Measurements and Theo- 
retical Calculations 

In Savin et al. (2002) we presented our results for DR measurements and calculations of 
oxygen-like Fe xix forming Fe xvm. Measurements were carried out using the heavy ion Test 
Storage Ring (TSR) at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. 
We also calculated the relevant DR resonance strengths and energies using our state-of-the- 
art AUTOSTRUCTURE code which uses a perturbative MCBP method (Badnell 1986). 
Overall reasonable agreement is found between experiment and our theoretical calculations. 
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Figure 1: Fe XIX DR laboratory rate coefficients due to 2 -4 2 (e~ 4- 2s 2 2p 4 — > 2£ 6 nl') and 
2 — > 3 (e~ 4 - 2s 2 2p 4 -4 2fi3lnl') core excitation; the latter are clearly seen in the collisionally 
ionized zone ( kT « 300 — 1000 eV, Mazzotta et al. 1998) as the 3131' ( kT w 400 eV), 314V 
( kT ~ 600 eV), 313V (kT ~ 750 eV), etc., manifolds: (a) the TSR experiment (Savin et al. 
1999, 2002), and (b) the present MCBP calculations. 

We have used our measurements to produce a Maxwellian-averaged rate coefficients for DR 
of Fe xix. Our experimentally-derived rate coefficient is estimated to be good to better than 
20%. Our AUTOSTRUCTURE results are in excellent agreement with our experimental 
rate coefficient, as seen in Fig. 1. However, these investigations turn up inaccuracies in the 
earlier recommended data by Mazzotta et al. (1998), Roszman (1987), and Jacobs et al. 
(1977), as seen in Fig. 2, indicating the utility of the newly computed data. 


2.4 Calculations of Fluoresence and Auger Yields 

In Gorczyca et al. (2003), we have investigated the accuracy of the fluoresence and Auger 
yields due to innershell ionization of boron-like ions and neon-like ions producing a Is hole, 
thereby creating, respectively, beryllium-like and fluorine-like ions which can then either 
radiatively stabilize or autoionize. We have carried out calculations up to and including zinc. 
The currently used Is fluoresence data base currently being used is based on 30 year old 
atomic calculations and simplistic isoelectronic extrapolations. Our results provide reliable 
new atomic data for two important isoelectronic sequences and cast doubt on the reliability 
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Figure 2: Fe XIX DR Maxwellian rate coefficients: (a) Thick solid curve - experimental data, 
dotted curve - MCBP calculations, thin solid curve - RR rate coefficient, (b) solid curve 
- present results, dot-dashed curve - Jacobs et al. (1977), dot-dot-dashed curve - Roszman 
(1987), dashed curve - recommended data of Mazzotta et al.(1998), dotted curve - MCDF 
(Savin et al. 2002). 

of the Is fluoresence database currently being used for the isoelectronic sequences of other 
second row elements. 


2.5 Calculated DR data for dynamic finite-density plasmas 

In Badnell et al. (2003) we outlined a program for the assembly of a comprehensive dielec- 
tronic recombination database within a generalized collisional-radiative (GCR) framework. 
It is valid for modeling ions of elements in dynamic finite-density plasmas. The resolution 
and precision of the data are tuned to spectral analysis and so are sufficient for prediction of 
the dielectronic recombination contributions to individual spectral line emissivities. The fun- 
damental data are structured according to the format prescriptions of the Atomic Data and 
Analysis Structure (ADAS) and the production of relevant GCR derived data for application 
is described and implemented following ADAS. The requirements on the dielectronic recom- 
bination database are reviewed and the new data are placed in context and evaluated with 
respect to older and more approximate treatments. We have benchmarked our results using 
laboratory measurements where available and found good agreement. All data produced are 
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Figure 3: DR cross sections for Ne 2+ for 2s — » 2p core excitation, (a) Comparison of cross 
sections obtained with R-matrix ( solid curve) and autostructure ( dotted curve) codes. 
The cross sections were convoluted with a Gaussian electron distribution of 0.1 eV FWFM. 
(b) Contribution of different decay modes to the total DR cross section. The region below 
21 eV is devoid of resonances formed due to a 2s — > 2p core excitation. 

publicly available via the website http://homepages.Hmich.edu/~gorczyca/drdata. 


2.6 Calculations of DR and RR Data for the Oxygen Isoelectronic 
Sequence 

In Zatsarinny et al. (2003), we have calculated DR and RR data for oxygen-like ions forming 
fluorine-like ions as part of the assembly of a level-resolved DR and RR database necessary 
for modeling of dynamic finite-density plasmas (see above). In this paper we presented total 
DR and RR rate coefficients for F + to Zn 22+ and discussed the results. By comparison 
between perturbative and R-matrix results, we find that RR/DR interference effects are 
negligible even for the lowest-charged (where interference is most likely) F + member (see 
Fig. 3). These RR and DR data are suitable for modeling of plasmas under conditions of 
collisional ionization equilibrium and non-equilibrium ionization. 






Figure 4: Maxwellian-average DR rate coefficients for Fe 20+ as a function of electron tem- 
perature (in eV). The thick solid curve represents our MCBP calculations, the solid squares 
represent the experimentally derived rate coefficients for just the 2 2 core excitations 

(Savin et al 2003). Also shown are our contributions due to just the 2 —> 2 core excitations 
( dotted line), the 2-4 3 core excitations (dash- dot- dotted line), and RR ( thin solid curve). 
Comparison is given with the LS-coupling calculations of Jacobs et al. (1977, solid circles), 
the unpublished LS-calculations of Rozmann as given by Arnaud & Raymond (1992, dot- 
dashed curve), and the recent calculations of Gu (2003a, short-dashed curve). For reference 
we have plotted up the temperature ranges where Fe 20+ is predicted to form in photoion- 
ization equilibrium (Kallman & Bautista 2001) and collisional ionization equilibrium (Gu 
2003b). 

2.7 Calculations of DR and RR Data for the Carbon Isoelectrbnic 
Sequence 

In Zatsarinny et al. (2004a), we have calculated DR and RR data for carbon-like ions forming 
nitrogen-like ions as part of the assembly of a level-resolved DR and RR database necessary 
for modeling of dynamic finite-density plasmas (see above). In this paper we presented total 
DR and RR rate coefficients for N + to Zn 24+ and discussed the results - see fig. 4 for a 
comparison between our theoretical results, the experimental ones of our collaborator Dr. 
Daniel Wolf Savin, and earlier theoretical results the accuracy of which were compromised 
by the computational limitations at that time and the resulting theoretical approximations 
required to make the calculations tractable. These RR and DR data are suitable for mod- 
eling of plasmas under conditions of collisional ionization equilibrium and non-equilibrium 
ionization. 
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Figure 5: Behavior of Maxwellian averaged Figure 6: Computed DR resonance struc- 
DR rate coefficients for the fluorine-like se- ture for neon-like ions Al 3+ , Ar 8+ , and 
quence as a function of atomic number. Fe 16+ . These data are subsequently con- 
Whereas the rate coefficients for 2 — ► 3 voluted analytically to yield Maxwellian- 
DR scale smoothly with Z (a), the rate averaged DR rate coefficients, which are 
coefficients for 2 — > 2 DR show irregular then condensed further into simple formu- 
dependence on Z (b). las using fitting coefficients. 

2.8 Calculations of DR and RR Data for the Fluorine and Neon 
Isoelectronic Sequences 

We also completed studies for the fluorine-like and neon-like sequences. These data have been 
made available on our web page (http://homepages.wmich.edu/~gorczyca/drdata). In 
Fig. 5, we show our computed results for the entire sequence (Zatsarinny et al., in prepa- 
ration). For Ne-like ions, where no experimental data exists, typical detailed resonance 
structure for selected ions is shown in Fig. 6 (Zatsarinny et al., 2004b). 
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Thank you for your continued support. 
Sincerely, 


Thomas W. Gorczyca and Kirk T. Korista 
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